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Field of the Invention 

The invention relates generally to the field of signal analysis, and more 
particularly, to a system and method for detecting the frequency, amplitude and/or phase 
of one or more tones comprised within an input signal. 

Description of the Related Art 

In many applications, it is necessary or desirable to precisely locate one or more 
tones comprised in a signal. This need arises in many fields including 
telecommunications, radar, sonar, networking, measurement, and various other 
applications. Prior art techniques for detecting tones in a signal generally may not 
produce accurate results and/or may suffer from slow performance. 

The discrete Fourier transform (DFT) is a popular tool for analyzing signals. 
However, before an input signal is transformed, it is quite often windowed with a 
windowing function. (It is noted that the action of capturing a finite-length sequence of 
samples of the input signal automatically implies a rectangular windowing.) The 
frequency transform F(n) of the windowed input signal will typically exhibit multiple 
scaled and shifted versions of transform function W, i.e., the transform of the window 
function. Each sinusoidal component of the input signal expresses itself as a pair of such 
shifted versions, one version shifted up to the frequency fj of the sinusoidal component, 
and the other shifted down to frequency -f y The positive frequency version is referred to 
herein as a positive frequency image, and the negative frequency version is referred to 
herein as a negative frequency image. When a sinusoidal component frequency fj is small 
compared to the sample rate, the positive frequency image and the negative frequency 
image for the sinusoidal component may overlap in frequency space. Similarly, when a 
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sinusoidal component frequency fj is close to one-half the sample rate, the positive 
frequency image and the negative frequency image for the sinusoidal component may 
overlap. Furthermore, when two sinusoidal components have frequencies that are close 
together, their positive images and negative images may overlap. 

Prior art techniques for tone estimation quite often focus on identifying the peaks 
in the magnitude spectrum |F(n)|. The peaks roughly determine the frequency of the 
corresponding tones. However, because of the cross-interaction of the images from other 
tones, or the negative frequency image from the same tone, the peak of a positive 
frequency image may be perturbed away from a purely scaled and frequency-shifted 
version of the template function W. Thus, parameter estimation techniques which 
compute parameters for a given tone based only on transform array values (i.e. DFT 
values) in the vicinity of a corresponding image peak may not produce accurate results. 
Therefore, there exists a substantial need for a system and method which could estimate 
tone parameters from the transform array with increased accuracy. 
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Summary of the Invention 

The present invention comprises various embodiments of a system and a method 
for estimating signal parameters (e.g., one or more of frequency, amplitude and/or phase) 
of one or more tones present in an input signal. More particularly, one embodiment of 
the invention comprises a system and method for estimating parameters for a tone based 
on a frequency transform F(n) of the input signal. The input signal may be windowed 
with a window function w(n) and transformed into the frequency domain. 

The tone in the input signal may express itself in the frequency domain as an 
additive combination of two spectra, one centered at the tone frequency and the other at 
the negative of the tone frequency. These two spectra are referred to herein as the 
positive frequency image and the negative frequency image respectively. A tone in the 
input signal may also be affected by spectra (e.g., positive or negative images) from other 
tones present in the signal. 

The method may comprise first receiving samples of the input signal, wherein the 
input signal includes the tone. A frequency transform of the samples may then be 
generated, e.g., a Fourier transform of the samples. The method may then identify a 
frequency location proximate to an amplitude peak in the frequency transform, wherein 
the amplitude peak corresponds to the tone. Two or more frequency bins may then be 
selected proximate to the identified frequency location in the frequency transform. The 
method preferably selects one or more frequency bins located on either side of the 
frequency location of the amplitude peak. In other words, the method may select one or 
more bins located on one side of the amplitude peak and one or more bins located on the 
other side of the amplitude peak. 

The method may then determine a tone frequency value that minimizes a 
difference between two or more expressions, e.g., at least a first expression and a second 
expression. Stated another way, the method may select a tone frequency value that makes 
a plurality of expressions most nearly equal. 
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Each of the plurality of expressions may comprise a sum of one or more 
numerator terms divided by a sum of one or more denominator terms, wherein: 

1) each of the plurality of expressions includes a tone frequency variable; 

2) each numerator term and each denominator term corresponds to one of the 
frequency bins; and 

3) a ratio of each numerator term and its corresponding denominator term 
represent a complex amplitude of the tone at a respective bin. 

Thus, for example, the plurality of expressions may comprise a first expression 
and a second expression having the above characteristics. 

In each expression, the tone frequency variable may represent a correct tone 
frequency value of the tone. The first expression may be approximately equivalent to the 
second expression when the correct tone frequency value is used for the tone frequency 
variable in the first and second complex expressions. 

The step of determining a tone frequency value may comprise computing a 
plurality of differences between the first expression and the second expression for 
different respective tone frequency values of the tone frequency variable, and then 
selecting the tone frequency value that produces a smallest difference. The method may 
compute the plurality of differences by performing a Newton-Rhapson root finding 
method. 

The expressions may be real expressions or complex expressions. Where the 
expressions are complex expressions, the method may involve minimizing a difference 
between an amplitude of the first complex expression and an amplitude of the second 
complex expression. 
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In the preferred embodiment, the first expression and the second expression have 
the form: 
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wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable fj. 

When the frequency transform of the samples comprises generating a power 
spectrum of the samples, the first expression and the second expression have the form: 
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wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable f^ 

The determined tone frequency value may comprise an approximation of the 
correct frequency of the tone. In other words, the determined tone frequency value 
represents an approximation of the correct tone frequency value. The method described 
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herein operates to find a more exact approximation of the frequency of the tone than prior 
art techniques, even in the presence of noise or interference from other tones, or the 
negative image of the respective tone being found. 

After the tone frequency value has been determined, the method may comprise 
storing the determined tone frequency value in a memory, or outputting the determined 
tone frequency value, e.g., on a display or to other systems for further analysis or 
processing. 

The method may also then compute one or more of the amplitude and phase of the 
tone using the determined tone frequency value. For example, once the ft value is known, 
the amplitude and phase can be computed as: 
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When the frequency transform of the samples comprises generating a power 
spectrum of the samples, the amplitude can be computed as: 
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Where the input signal comprises a plurality of tones, the method may operate to 
correctly find signal parameters of any one of the plurality of tones. Where it is desired 
to locate two or more of, or all of, the tones present in the input signal, the above steps 
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may be iteratively performed for each of the plurality of tones to determine at least one 
signal parameter for each of the plurality of tones in the input signal. The above method 
may be performed iteratively wherein, after a respective tone is identified, the respective 
tone may be subtracted from the input signal, and the above process repeated for another 
5 tone, and so on. 

Thus the method of the present invention operates to find more correct 
approximations of signal parameters of tones present in an input signal, and operates 
more efficiently than prior art techniques. 

10 
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Brief Description of the Drawings 

A better understanding of the present invention can be obtained when the 
following detailed description of the preferred embodiment is considered in conjunction 
5 with the following drawings, in which: 

Figure 1A illustrates a system configuration 100 for determining the signal 
parameters associated with one or more sinusoidal tones comprised within an input 
signal; 

Figure IB illustrates one embodiment for tone detection system 120; 
10 Figure 2A illustrates one embodiment of tone detection system 120 comprising a 

computer-based measurement or data acquisition system, where signals generated by 
signal reception device SRD are presented to computer 102 through signal conditioning 
system 108 and data acquisition (DAQ) device 104; 

Figure 2B illustrates a second embodiment of tone detection system 120 
15 comprising a computer-based measurement system, where signals generated by signal 
reception device SRD are presented to computer system 102 through data acquisition 
(DAQ) device 104; 

Figure 3 presents a flowchart for one embodiment of a tone detection system 
according to the present invention; 
20 Figure 4 illustrates a windowing operation being performed on an input signal to 

generated a windowed input signal; 

Figure 5 illustrates the magnitude of transform array F(k) for a typical windowed 
input signal comprising a single sinusoidal tone; 

Figure 6 illustrates a blowup of a generic magnitude peak 301 from the magnitude 
25 spectrum of Figure 5; and 

Figures 7A and 7B, 8A and 8B, 9A and 9B, 10A and 10B, and 11A and 11B 
illustrate front panels and block diagrams of a portion of a graphical program that 
implements one embodiment of the invention. 
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While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof are shown by way of example in the drawings and will 
herein be described in detail. It should be understood, however, that the drawings and 
detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the intention is to cover all modifications, equivalents and 
alternatives falling within the spirit and scope of the present invention as defined by the 
appended claims. 
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Detailed Description of the Embodiments 



Incorporation by Reference 

U.S. Patent Application Serial No. 09/753,164 titled "System and Method for 
5 Estimating Tones in an Input Signal", Filed December 27, 2000, is hereby incorporated 
by reference in its entirety as though fully and completely set forth herein. 

Figure 1A 

Figure 1A illustrates a system configuration 100 for performing signal processing 

10 on a signal comprising one or more tones. System configuration 100 may comprise a 
signal reception device SRD and a tone detection system 120. 

The SRD may be coupled to receive a signal from a device, unit under test (UUT) 
or a transmission medium 110, or any other system capable of transmitting a signal that 
may contain tones. The term "transmission medium" is used herein to refer generally to a 

15 device, unit under test (UUT) or a transmission medium 110 that may generate a signal 
including one or more tones. As used herein, the term "tone" includes a signal at a 
frequency, e.g., at a primary or single frequency, which may be contained within another 
signal or may itself be the signal. For example, the signal may simply comprise the tone 
and a small amount of noise. Alternatively, the signal may comprise the tone and one or 

20 more other signals, e.g., one or more other tones. As an example, a tone may comprise a 
sinusoidal signal having arbitrary amplitude, frequency and phase. 

As shown in Figure 1A, SRD may be coupled to a transmission medium 110. 
Transmission medium 110 may represent any of a variety of transmission media such as 
the atmosphere, free space, an optical fiber or fiber bundle, a communication bus (e.g. a 

25 network bus), a body of water or any other fluid, the earth, etc. In one embodiment, 
transmission medium 110 is the atmosphere, and signal reception device SRD comprises 
an antenna and a radio receiver. In a second embodiment, transmission medium 1 10 is a 
network bus connecting two or more computers, and signal reception device SRD is a 
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network interface card/board. In a third embodiment, transmission medium 110 is an 
optical fiber, and signal reception device SRD comprises an optical sensor. As noted 
above, element 110 may be any of various devices or mediums for generating or 
transmitting a signal 

5 Signal reception device SRD receives an input signal from the transmission 

medium or device 110 and may convert the input signal into a form suitable for 
presentation to tone detection system 120. The input signal may be electrical or non- 
electrical in nature. Signal reception device SRD may include analog-to-digital 
conversion hardware to digitize the input signal. Alternatively, analog-to-digital 

10 conversion hardware may be comprised within tone detection system 120. 

In one embodiment, signal reception device SRD may comprise a measurement 
device such as a microphone, an accelerometer, a spatial displacement sensor, a strain 
gauge, a pressure sensor, a temperature sensor (e.g., a thermocouple), a radiation sensor, 
an optical sensor, etc, or any combination thereof. In another embodiment, signal 

15 reception device SRD may represent an array of transducers or measurement devices of 
one or more types. SRD may thus be any of various transducers or sensors for receiving 
a signal. 

Tone detection system 120 may couple to signal reception device SRD to receive 
the input signal. In an alternate embodiment, the tone detection system 120 may be 

20 operable to receive the input signal from a memory, wherein the input signal was 
previously acquired and stored in the memory. 

Tone detection system 120 may be configured for detecting the frequency, 
amplitude and/or phase of one or more tones in the input signal. Tone detection system 
120 may comprise a processor or central processing unit (CPU) 140, memory 146, user 

25 input device(s) UID and a display device DD as shown in Figure IB. CPU 140 may be 
realized by any of a variety of computational devices such as a general purpose processor, 
a digital signal processor, a parallel processor, dedicated digital and/or analog circuitry, 
programmable gate array logic (e.g., an FPGA), etc., or any combination thereof. 
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Memory 146 may comprise any of a variety of memory devices such as random access 
memory (RAM) and/or read-only memory (ROM), as described further below. Tone 
detection system 120 may also include specialized data acquisition and/or signal 
conditioning hardware, interface hardware, etc., or any combination thereof. 
5 Tone detection system 120 may comprise any of various devices, such as a 

programmable computer system, a computer-based system such as a VXI-based system, a 
PXI-based system, a GPIB-based system, a computer-based data acquisition system, or a 
dedicated test instrument, such as a dynamic signal analyzer, an oscilloscope or any other 
signal acquisition and/or analysis device. 

10 Tone detection system 120 may operate on samples of the input signal X 

generated by signal reception device SRD and thus may identify the frequency, phase 
and/or amplitude of one or more tones in the input signal. The frequency, phase and/or 
amplitude of the one or more tones may be presented to a user through the display device 
DD or some other output device, may be further processed by a software program 

15 executing in the tone detection system 120, may be provided to another system for 
processing, and/or may be stored to memory for future use. 

User input device(s) UID may comprise a keyboard, a pointing device such as a 
mouse or trackball, a touch pad (such as those used in modern laptop computers for 
cursor control), a touch sensitive display screen, etc., or other input devices. In one 

20 embodiment, user input device(s) UID may include use of a graphical user interface or 
panel configured with various control icons such as buttons, knobs, sliders, switches, 
indicators, etc., or any combination thereof. A user provides input to tone detection 
system 120 through user input device(s). Tone detection system 120 may manage a 
graphical user interface through display device DD and user input device(s) UID. 

25 

Figures 2A and 2B 

Figure 2A and 2B illustrate exemplary embodiments of tone detection system 
120. As shown, tone detection system 120 may comprise a computer 102, a data 
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acquisition (DAQ) device 104 coupled to the computer 102, and optionally a signal 
conditioning system 108 coupled to the DAQ device 104. Signal reception device SRD 
may comprise transducers, sensors, and/or receiving devices that couple to DAQ device 104 
through the signal conditioning circuitry 108. 
5 As shown, signal reception device SRD is configured and/or coupled to acquire 

signals from the transmission medium 110. The input signals acquired by signal reception 
device SRD may be optionally conditioned by the signal conditioning system 108 as shown 
in Figure 2 A. The conditioned input signals may then be provided to DAQ device 104 as 
shown. Signal conditioning system 108 may connect to DAQ device 104 via one or more 
10 cables. 

In one embodiment, signal conditioning system 108 may comprise an external 
chassis 122 housing one or more signal conditioning modules 124 and optionally terminal 
blocks 126. Signal conditioning system 108 may be used to perform signal conditioning on 
the signals generated by signal reception device SRD. As used herein, the term "signal 

15 conditioning" may include one or more of amplifying, linearizing, limiting, isolating, 
filtering, switching and/or multiplexing signals), among other signal processing functions. 
Signal conditioning system 108 may advantageously reduce the level of noise in the signals 
transmitted to DAQ device 104. DAQ device 104 may receive conditioned signals from 
signal conditioning system 108 as shown in Figure 2A. Alternatively, DAQ device 104 

20 may directly receive the input signal from signal reception device SRD as shown in Figure 
2B. DAQ device 104 may operate to perform analog to digital (A/D) conversion and 
provide the resultant digital signals to computer 102 for processing. 

Computer system 102 may include various standard components, including a 
processor or central processing unit (CPU) 140, system memory 146, non-volatile 

25 memory, one or more buses, and a power supply. DAQ device 104 may be a specialized 
system for acquiring digital and/or analog signals from external devices. Thus, DAQ 
device 104 may include analog to digital (A/D) conversion circuitry and/or digital to 
analog (D/A) conversion circuitry. Examples of the DAQ device 104 include "E series" 
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DAQ boards from National Instruments Corporation. DAQ device 104 may also 
comprise a computer-based instrument board, such as an oscilloscope, a digital multimeter 
(DMM), a dynamic signal analyzer, an arbitrary waveform generator, etc. 

In one embodiment, computer 102 may comprise input/output (I/O) slots into 
5 which DAQ device 104 may be coupled. In another embodiment, computer 102 may 
comprise a VXI (VME Extensions for Instrumentation) chassis and bus, a GPIB (General 
Purpose Interface Bus) interface card, a serial port or parallel port by which DAQ device 
104 may be coupled to the computer 102. 

Tone detection system 120, e.g., computer system 102, preferably includes at least 

10 one memory medium on which computer programs according to the present invention may 
be stored. The term "memory medium" is intended to include various types of memory or 
storage, including an installation medium, e.g., a CD-ROM, or floppy disks 104; a 
computer system memory or random access memory such as DRAM, SRAM, EDO RAM, 
Rambus RAM, EPROM, EEPROM etc.; or a non-volatile memory such as a magnetic 

15 media, e.g., a hard drive, or optical storage. The memory medium may comprise other 
types of memory as well, or combinations thereof In addition, the memory medium may 
be located in a first computer in which the programs are executed, or may be located in a 
second different computer which connects to the first computer over a network. In the latter 
instance, the second computer may provide the program instructions to the first computer 

20 for execution. Also, the computer system 102 may take various forms, including a personal 
computer system, mainframe computer system, workstation, network appliance, Internet 
appliance, personal digital assistant (PDA), television system, dedicated test or 
measurement instrument or other device. In general, the term "computer system" can be 
broadly defined to encompass any system having a processor which executes instructions 

25 from a memory medium. 

The memory medium preferably stores a software program according to one 
embodiment of the present invention for detecting one or more tones in the input signal. 
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More particularly, the software program may be operable to analyze the input signal to 
determine the frequency, phase and amplitude of one or more tones in the input signal. 

The software program may be implemented in any of various ways, including 
procedure-based techniques, component-based techniques, object-oriented techniques, or 
5 neural net based learning techniques, among others. For example, the software program 
may be implemented using ActiveX controls, C++ objects, Java objects, Microsoft 
Foundation Classes (MFC), or other technologies or methodologies, as desired. A 
processor, such as the host CPU, executing code and data from the memory medium, or a 
programmable device configured according to a net list, may comprise embodiments of a 

10 means for determining the frequency, phase and amplitude of the one or more tones 
embedded in the input signal according to the methods described below. 

As noted above, the tone detection system 120 may also or instead include 
reconfigurable logic or programmable logic, such as an FPGA, that implements the 
methods described herein. 

15 Various embodiments further include receiving, storing, and/or transmitting 

instructions and/or data implemented according to the present invention upon a carrier 
medium. Suitable carrier media include a memory medium as described above, as well as 
signals such as electrical, electromagnetic, or digital signals, conveyed via a 
communication medium such as networks and/or a wireless link. 

20 

Figure 3 - Determining Signal Parameters of a Tone 

Figure 3 illustrates one embodiment of a method for determining one or more 
signal parameters, e.g., one or more of the frequency, amplitude and/or phase, of one or 
more tones present in the input signal. In one aspect, the method may be referred to as an 
25 aliasing compensation method for compensating for the effects of the negative image of a 
tone on the positive image of the tone. In another aspect, the method may be used to 
more correctly identify signal parameters of one or more tones in a signal, compensating 
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for the effects of the negative image of the tone, and the positive and/or negative images 
of other tones present in the signal. 

The method of Figure 3 may be implemented by execution of a computer program 
stored on the memory medium as described above. The method may also be 
5 implemented in programmable logic, such as an FPGA. The method may also be 
implemented by a combination of hardware and software, as desired. It is noted that 
various steps may be performed concurrently or in a different order than that shown, 
and/or some steps may possibly be omitted, as desired. 

In step 202 the method may receive samples x(n) of the input signal, wherein the 
10 input signal includes the tone (or one or more tones). The samples x(n) may be provided 
by signal reception device SRD or may be received from a memory medium, e.g., having 
been previously recorded/captured from signal reception device SRD. Alternatively, the 
input signal samples may be simulated samples generated by a simulator (e.g. a CPU 
executing simulation code). Various embodiments of the invention contemplate a wide 
15 variety of possible sources for the input signal samples x(n). 

The input signal may comprise a single sinusoidal tone in the presence of noise. 
Thus, the input signal and/or the tone may be modeled by the expression: 
x{ri) = A t oos(27f i n + <p t ) 
Where Af is the amplitude of the tone,// is the normalized frequency of the tone, 
20 n is from 0 to N-l, and (pi is the initial phase of the tone. 

In step 204 the CPU 140 may multiply the input samples by a known window 
function w(n) to generate a windowed input signal y(n) = w(n) * x(n) as suggested by 
Figure 4. Alternatively, the CPU 140 may receive, e.g., in step 202, a discrete set of 
25 samples, which effectively form a windowed version of the signal. The window function 
w(n) may have any of a variety of forms. For example, the window function may be a 
rectangular window, a triangular window, a raised cosine window, a Hanning window, 
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etc. The window function may, for example, be any type of cosine window, such as a 
Harming, Hamming, Blackman Harris, and Flat Top window, among others. 

In step 204, the method may also perform a frequency transform, such as the 
discrete Fourier transform (DFT), e.g., an FFT (Fast Fourier Transform), on the 
5 windowed input signal y(n) to generate a transform array F(n), where n is a frequency bin 
index which may range from 0 to N-l, or any interval of length N, where N is a positive 
integer. In one embodiment or one example, the transform array F(n) may be modeled by 
the transform of the sinusoidal tone, i.e., 

F(n) = A?'* W(f n + /, ) + V* w (fn ~ fi ) 
10 where W(f) represents the Fourier transform of the window w(n). It is noted that the 
relationship between frequency f and frequency bin number k is given by 

f=f s *(k/N), 

where f s is the sample rate. The magnitude of the window transform W(f) typically has 
even symmetry and attains a maximum at f = 0. Thus, the function W(f-Q attains a 

15 maximum magnitude at frequency f = f 1? and the function W(f+f { ) attains a maximum 
magnitude at frequency f = -f i( The first term in the expression above, i.e., 

P(f) = (A/2)expO»W(f-f i ), 
is referred to herein as the "positive-frequency image" since its center frequency occurs at 
the positive frequency f { . The second term in the expression above, i.e., 

20 N(f) - (A/2)exp(-j <p )W(f+Q 

is referred to herein as the "negative-frequency image" since its center frequency occurs 
at the negative frequency -f { . Thus, the transform array F(n) may include a positive- 
frequency image and negative-frequency image which combine additively (in the sense of 
complex addition). The input signal may also include noise, other spurious tones, or 

25 other valid tones. 

If tone frequency f { stays away from zero or f s /2, and/or the sample size N is 
sufficiently large, the overlap between the positive and negative frequency images may be 
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small, and thus, their individual identities may be apparent in the transform array F(n). 
The magnitude function |F(n)| will thus exhibit two peaks which correspond to the 
positive and negative frequency images. The frequency locations of one of these peaks 
(i.e., the peak that occurs in the range of positive frequencies) may be used as an estimate 
5 for the tone frequency fj. 

Conversely, if the tone frequency is close to zero or f s /2, and/or, the sample size N 
is sufficiently small, the positive-frequency image and negative frequency image may 
overlap significantly. Thus, their individual identities may not be apparent in the 
transform array F(n). In other words, transform array F(n) restricted to positive 
10 frequencies may be a poor approximation to the positive frequency image. Thus, the 
frequency location at which the magnitude function |F(n)| attains a maximum, when 
considered over positive frequencies, may only be a crude initial approximation to the 
tone frequency fj. 

Figure 5 is a plot of the magnitude of transform F(n) corresponding to a typical 
15 windowed input signal x(n). Note that the transform F(n) has a symmetry given by F(n) 
= F(k+N) for any integer k. In particular, F(-k)=F(N-k). Thus, frequency bin numbers 
between N/2 and N may be interpreted as negative frequencies. 

In step 206, method may identify a frequency location proximate to an amplitude 
peak in the frequency transform, wherein the amplitude peak may correspond to the tone. 
20 For example, in step 206 the method may determine the maximum amplitude peak in the 
input signal, which is presumed to be the peak of the tone. In one embodiment, the 
method may scan the DFT magnitude values |F(n)| over the range of positive frequency 
bins to determine the potential local positive maximum amplitude, e.g., to determine the 
bin index k which achieves the maximum magnitude. Where multiple tones are present 
25 in the input signal, the method may find multiple amplitude peaks and select one for 
subsequent processing. 

In step 208 the method may select two or more frequency bins proximate to the 
identified frequency location in the frequency transform. Thus the method may select 
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two or more bins proximate to (or around) the positive maximum determined in 206. The 
method preferably selects frequency bins located on either side of the frequency location 
of the amplitude peak. In other words, at least one of the two or more frequency bins is 
on each side of the frequency location, i.e., one or more bins on a first side of the 
5 frequency location of the amplitude peak and one or more bins on the other side of the 
frequency location of the amplitude peak. In one embodiment, in step 208 the method 
may select several bins k that are in the neighborhood of k,^, wherein k max is the integer 
bin index value k in the range from 0 to N/2 which maximizes the magnitude of F(n). 

It is noted that the maximum of magnitude function |F(n)| considered as a function 

10 of continuous frequency typically does not occur at the integer value k^, although it 
should occur somewhere in the interval between and k 2 Figure 6 illustrates a blowup 
of the positive frequency magnitude peak 301 in the neighborhood of bin index value 
k^. In the example of Figure 6, the method may select - 1 and k^ + 1 as the bins 
proximate to the frequency producing the positive maximum amplitude. In the preferred 

15 embodiment, the method selects three frequency bins for processing. However, the 
method may select 2, 3, 4, or a greater number of bins, as desired. 

In step 210 the method may determine a tone frequency value that minimizes a 
difference between two or more expressions, e.g., at least a first expression and a second 

20 expression. Stated another way, the method may involve selecting a tone frequency value 
that makes a plurality of expressions most nearly equal. 

Each of the plurality of expressions may comprise a sum of one or more 
numerator terms divided by a sum of one or more denominator terms, wherein each of the 
plurality of expressions includes a tone frequency variable, wherein each numerator term 

25 and each denominator term corresponds to one of the frequency bins, and wherein a ratio 
of each numerator term and its corresponding denominator term represent a complex 
amplitude of the tone at a respective bin. 
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Thus, for example, the first expression may comprise a sum of two or more 
numerator terms divided by a sum of two or more denominator terms, wherein the first 
expression includes a tone frequency variable. In a similar manner, the second expression 
may comprise a sum of one or more of the numerator terms divided by a sum of one or 
5 more of the denominator terms, wherein the second expression includes the tone 
frequency variable. The number of numerator terms and corresponding denominator 
terms in the first expression are preferably different than the number of numerator terms 
and corresponding denominator terms in the second expression. 

In each expression, each numerator term and its corresponding denominator term 

10 correspond to one of the frequency bins. In one embodiment, a ratio of each numerator 
term and its corresponding denominator term represent a complex amplitude of the tone 
at a respective bin. In each expression, the tone frequency variable may represent a 
correct tone frequency value of the tone. The first expression may be approximately 
equivalent to the second expression when the correct tone frequency value is used for the 

15 tone frequency variable in the first and second expressions. 

As one example of the expressions used, where two frequency bins are selected, 
the first expression may comprise a sum of two numerator terms divided by a sum of the 
corresponding two denominator terms, wherein the two numerator terms and the two 
denominator terms correspond to the two frequency bins. The second expression may 

20 comprise one numerator term divided by one corresponding denominator term 
corresponding to one of the two frequency bins. 

As another example, where three frequency bins are selected, the first expression 
may comprise a sum of three numerator terms divided by a sum of the corresponding 
three denominator terms, wherein the three numerator terms and the three denominator 

25 terms correspond to the three frequency bins. The second expression may comprise a 
sum of two of the three numerator terms divided by a sum of two of the three 
corresponding denominator terms, corresponding to two of the three frequency bins. 
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The step of determining a tone frequency value may comprise computing a 
plurality of differences between the first expression and the second expression for 
different respective tone frequency values of the tone frequency variable, and then 
selecting the tone frequency value that produces a smallest difference. The method may 
5 compute the plurality of differences by performing a Newton-Rhapson root finding 
method, as is known in the art. 

The expressions may be real expressions or complex expressions. Where the 
expressions are complex expressions, the method may involve minimizing a difference 
between an amplitude of the first complex expression and an amplitude of the second 
10 complex expression. 



In one embodiment, the first expression and the second expression have the form: 



15 



20 



k+M , x 

Y\F{n)W\f n -f l )-F'(n)W{f u + /,)) 

Jfe+1 



k+M/ \ 

U w (f»-fit-\ w (f» + fif) 

k+l 



k+M-l v 

Y\F(n)W\f n -fd-F\n)W{j ,+/,)) 

k+l 



k+M-l / 

k+l 



wherein: 



F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable fj. 



Introducing the function G(f): 
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G(f) = 



k+M, A 
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k+l 
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k+M-l 
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+M-1 k+M/ \ 
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This results in: 
G(f) = 0 

5 This equation is a function of f, and this equation can be solved by applying 

Newton Raphson ? s root finding algorithm to determine the tone frequency^*. The theory 
behind the use of these expressions is discussed further below. 

When the frequency transform of the samples computed in 204 comprises 
10 generating a power spectrum of the samples, the first expression and the second 
expression have the form: 



k+M J 



k+l 



k+M 



5y(/,-/,)f 



k+l 



k+M-l 
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k+l 
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k+i 



15 wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable f { . 



20 



Let: 
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k+M k+M-l k+M-l k+M 

G(/)=Z^(")|X|^-/)|)X I Wn-ff- Un n P\ W (L-f)\)xB W (fn-ff 

k+l k+l k+l k+l 

This results in: 



G(f) = 0 

5 Again, Newton Raphson root finding may be performed to solve for the frequency 

of the tone. The theory behind the use of these expressions is discussed further below. 

The tone frequency value determined in step 210 may comprise the correct 
frequency of the tone. As used herein, the term "correct" includes correct or substantially 

10 correct. In other words, the determined tone frequency value represents a good 
approximation of the correct tone frequency value. The method described herein operates 
to find a more exact approximation of the frequency of the tone, even in the presence of 
noise or interference from other tones (i.e., the positive and/or negative images from other 
tones), or the negative image of the respective tone being found. The method is also 

1 5 more computationally efficient than prior art techniques. 



Where the input signal comprises a plurality of tones, the method may operate to 
correctly find signal parameters of a first tone of the plurality of tones. Where it is 
desired to locate two or more of, or all of, the tones present in the input signal, steps 206 
20 - 210 may be iteratively performed for each of the plurality of tones to determine at least 
one signal parameter for each of the plurality of tones in the input signal. 

After the tone frequency value has been determined in step 210, in step 212 the 
method may then optionally compute one or more of the amplitude and phase of the tone 
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using the determined tone frequency value. For example, wherein the FFT spectrum is 
available, once the fj value is known, the amplitude and phase can be computed as: 



k+M, , 

Y\F{nW\f n -fi)-F\n)W{f n +f i )) 



Afto =C = 



k+M ( ? i \ 

Zrc/,-/,)| f-m/n+M) 



k+i 



Where the FFT power spectrum is available, once the frequency is known, the 
amplitude can be computed as: 



k+M 

^(n)\x\W(f n -/ ( )|) 



4 = 



k+\ 




k+l 



In step 214 the method may then comprise storing the determined tone frequency 
value, as well as the amplitude and/or phase if these are computed, in a memory, or 
outputting the determined tone frequency value, amplitude and/or phase, e.g., on a 
display. 



Theory of Operation of the Expressions 

The following describes the derivation and theory behind the expressions used in 
the preferred embodiment. 

As noted above, given a single tone signal in the time domain, the input signal 
20 comprising the tone can be expressed in following form: 



Where Aj is the amplitude,/^ is the normalized frequency, n is from 0 to N-l and 
q>i is the initial phase. Here the subscript "i" is used to represent a single tone in what 



x(n) = 4- cos(2^/j.w + q> x ) 



(1) 
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could be a multitone signal. 



If the FFT is computed based on the windowed version of x(n) as performed in 
step 204, then in the frequency domain, as noted above, each FFT bin can be computed as 
5 follows: 

F(/i) = +f i ) + A i e j «W(f n -f t ) (2) 



In the above equation, the first term of the equation represents the negative image 
of the tone and the second term represents the positive image of the tone. 
10 Each equation can then be converted to a different form. For example, given the 

following: 

C = V* (3) 

Then 

F(n) = C*W(f n + /,) + CW(f n - /,) (4) 



15 



Here, F(n) is the wth value of the single sided scaled FFT spectrum. 
The windowed frequency response can be easily computed numerically once the 
window coefficients are known. 



20 Estimate Amplitude/Frequency/Phase based on FFT spectrum 

If the FFT spectrum is available, each particular single tone can be roughly 
determined by finding the local maximum in the sense of amplitude. Here assume 
F(k+1), F(k+2) v . F(k+M) are the M bins around the single tone. This results in a set of 
nonlinear equations of : 

25 

Fin) = C*W(f n + /J) + CW(f n -f t ) (5) 
Where n - k+1, k+2, . . . k+M; 
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From this point, the following equation can be derived: 
_ F(n)W*(f n -f i )-F\n)W(f n+ f i ) 
frif.-ftf-frif.+f.H 2 
For each n = k+1, . , .k+M. 

Furthermore, to reduce the numerical fluctuation, the equation (7) may be 
written as follows: 



k+M 

n F (nW'(f n -fd-F\n)W{f n +/,)) 

k+1 
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k+l 
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k+1 



+M-1 f \ 

I wVn-fii-^iL+ftf) 



This equation is a function of/J , and this equation may be solved in step 210 by 
applying Newton Raphson's root finding algorithm. 

Define G(f): 



G(f) = 



k+M , I k+M-l ( \ 

£(F(n)W\f n -f)-F*(n)W(f n+ f)]x £ W(f n - ftf -\W(f n + ff) 

k+\ | k+1 



(8) 
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YXnnW" (f n ~f)-F\n)W{f n +/)]x 

k+1 



k+M I \ 
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Now by solving the equation: 

G(f) = 0 (9) 
using Newton Raphson's root finding algorithm (or other suitable root find 
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methods), the tone frequency f{ can be determined 

Once thcff value is known, the amplitude and phase can be computed as follows: 



k+M y 



YjF{n)W\f n - fi )-F\n)W{f n +/,)) 
5 V* =C = ^Ss7 r— (10) 

un/n-f^-lnL+M 2 ) 



k+1 



Estimate Amplitude/Frequency based on FFT power spectrum 

If the FFT power spectrum is available, then the amplitude and frequency of the 
10 tone may be estimated based on the FFT power spectrum. Since the phase information is 
not available in the power spectrum, the effect of the negative frequency cannot be taken 
into account. The equation for computing the frequency are as follows: 

k+M 

2>(«)|x|r(/„-./;.)|) 

k+M 

k+1 

15 

Let equation (12) be: 

k+M k+M-1 k+M -I k+M 

G(f)=Z(\F(n)\x\W(f„-f)\)x X \W(f n -f)\ ~ S^(")|x|^(A-/)|) x ZI^-/)r 

k+1 k+1 k+l k+1 

This results in equation (13): 
20 G(/) = 0 

Again, the Newton Raphson root finding method may be applied to solve for the 
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(11) 
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frequency of the tone. 

Once the frequency is known, the amplitude can be computed by equation (13) as: 



4 = 



k+M 

U F (n)\x\W(f 



k+l 



k+M 
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k+l 



(14) 



In one embodiment, the power spectrum may be averaged before applying this 
tone estimation. The measurement is more stable when using the averaged power 
spectrum. 



10 Determining a Plurality of Tones 

The above method may be used to find a plurality of different tones present within 
a signal. For example, the above method may be applied iteratively to find each 
respective tone within the input signal. In one embodiment, as each tone is determined or 
found, the determined tone may be subtracted out from the input signal to produce a 

15 modified input signal, and the method may be then performed on the modified input 
signal to locate the next tone, and so on. 



Applications 

Embodiments of the present invention may be used in various applications. In 
20 general, embodiments of the present invention may be used in any system where it is 
desired to detect tones, e.g., sinusoidal tones, present in a signal, e.g., where it is desired 
to detect the precise frequency, amplitude and/or phase of the tones present in the signal. 
For example, an embodiment of the present invention may be used in a DTMF (Dual 
Tone Multi-Frequency) system for detecting tones present in a signal, such as a signal 
25 generated by a keypad of a telephone. Embodiments of the present invention are also 
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contemplated for use in applications involving sonar, radar (e.g. Doppler radar), 
frequency-shift keying applications, mechanical systems analysis, etc. For example, the 
reflections generated by multiple moving objects in response to a radar pulse have 
distinct frequencies dependent on their radial velocities with respect to the radar station. 
Thus, the frequencies of the reflections are usable for tracking the multiple moving 
objects. In another example, a mechanical system excited with a physical stimulus (e.g. 
an impulse) may manifest vibrations at one or more frequencies. The frequency, 
amplitude and/or phase of these vibrations may provide information to a system analyst 
about the nature of flaws in the mechanical system. Embodiments of the present 
invention may be used in a wide variety of applications, i.e. in any application where it is 
desirable to identify one or more tones present in an input signal. The above-mentioned 
applications are merely representative examples. 

Example Graphical Program 

Figures 7A -1 IB comprise screen shots of a portion of a graphical program which 
implements one embodiment of the invention. Figures 7A -1 IB comprise screen shots of 
a portion of a graphical program or Virtual Instrument (VI) written in the Lab VIEW 
graphical programming language. The graphical program is hierarchical in nature, where 
Figures 7 A and 7B comprise the top level front panel and block diagram, and figures 
8A/B, 9A/B, 10A/B, and 11A/B are various sub programs or sub-Vis of the program. 
Various other sub programs or sub-Vis of this graphical program that are not necessary to 
an understanding of this embodiment of the invention are not included for convenience. 

Figures 7A and 7B illustrate a top level front panel and block diagram, 
respectively, of the graphical program titled "estimate single tone". This graphical 
program includes a block titled "single tone info" ("single tone information"). 

Figures 8A and 8B illustrate a front panel and block diagram, respectively, of the 
graphical program block titled "single tone information" contained in the diagram of 
Figure 7B. This graphical program includes a block titled "tone est " ("tone estimation"). 
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Figures 9A and 9B illustrate a front panel and block diagram, respectively, of the 
graphical program block titled "tone estimation" contained in the diagram of Figure 8B. 
This graphical program includes two blocks labeled "G(f)" and "G'(f)" 

Figures 10A and 10B illustrate a front panel and block diagram, respectively, of 
5 the graphical program block titled "G(f)'\ 

Figures 11A and 11B illustrate a front panel and block diagram, respectively, of 
the graphical program block titled "G'(f)'\ 

Although the system and method of the present invention is described in 
10 connection with several embodiments, it is not intended to be limited to the specific 
forms set forth herein, but on the contrary, it is intended to cover such alternatives, 
modifications, and equivalents, as can be reasonably included within the spirit and scope 
of the invention as defined by the appended claims. 
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